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ABSTRACT 

We analyze high signal-to-noise spectra for 14 globular clusters (GCs) in the Som- 
brero galaxy, NGC 4594 = M104, obtained with the LRIS spectrograph on the Keck 
I telescope. We find generally good agreement between spectroscopic metailicities and 
previous photometric metallicites based on V —I colors. Further, we use measurements 
of various Balmer line indices on co-added spectra of metal-poor and metal-rich GCs to 
estimate ages for the GC subpopulations and find ages of 10-15 Gyr, with preferred ages 
around 11-12 Gyr and no detectable differences between the metal-rich and metal-poor 
subpopulations. However, because of model and measurement uncertainties, age differ- 
ences of a few (2-3) Gyrs cannot be ruled out. We also measure alpha to iron abundance 
ratios from Mg and TiO features and find enhancements on the order of [a/Fe] ~ +0.4, 
similar to the [a/Fe] abundance ratios observed in other old stellar populations. Finally, 
we combine our sample with radial velocities for 34 GCs from Bridges et al. (1997) and 
obtain virial and projected masses for the Sombrero of Myx = (8.1 ± 1-8) x 10 11 M Q 
and Mp = (5.3 ± 1.0) x 10 11 M Q within 17 kpc, respectively. Using 12 clusters within 
the central 4.5 kpc the corresponding numbers are Mvt = (2-1 ± 1.1) x 10 11 M Q and 
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M P = (2.0 ± 1.3) x 10 11 M G 

Subject headings: galaxies: star clusters — galaxies: abundances - galaxies: halos 
galaxies: bulges — galaxies: individual (NGC 4594) 



Based on data obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and the National Aeronautics and Space Administra- 
tion. 
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1. Introduction 

The Sombrero galaxy (M104 = NGC 4594) is well known for its dominating bulge/halo com- 
ponent which contributes with about 80% of the total luminosity (Kent 1988; Baggett, Baggett 
& Anderson 1998). This makes the Sombrero an important intermediary case between later- type 
spirals and elliptical/SO galaxies, and an ideal place to study the relation between spiral bulges and 
early-type galaxies. One interesting aspect of such studies is to compare the globular cluster (GC) 
systems of early- and late-type galaxies. The Sombrero has a very rich GC system with an esti- 
mated total of ~ 2000 clusters (Harris et al. 1984; Bridges & Hanes 1992; Forbes, Grillmair & Smith 
1997), about an order of magnitude more than the Milky Way and M31 (Harris 1996; Fusi Pecci 
et al. 1993), but comparable to many early- type galaxies of the same luminosity. Using archive 
WFPC2 images, Larsen, Forbes k, Brodie (2001, hereafter Paper I) detected a clearly bimodal 
color distribution for the Sombrero GCs, confirming hints by Forbes, Grillmair & Smith (1997). 
The number of red (metal-rich) and blue GCs were found to be roughly equal, in contrast to M31 
and the Milky Way where the metal-poor GCs constitute about 2/3 - 3/4 of the total population, 
but again similar to the situation in early-type galaxies. Thus, we suggested in Paper I that the 
metal-rich GCs in the Sombrero are associated with the bulge rather than with its disk, in which 
case a much more modest number of metal-rich GCs would have been expected. The WFPC2 data 
presented in Paper I also showed that the mean colors of the two GC populations in the Sombrero 
are very similar to those in early-type galaxies as well as the Milky Way, M31 and other spiral 
galaxies (Forbes, Brodie & Larsen 2001), with inferred mean metallicities of [Fe/H]~ —1.40 and 
-0.54. 

Although there may not be a one-to-one correspondence between star formation and formation 
of globular clusters, information about the age distribution of GCs would presumably provide valu- 
able insight about major star formation episodes in their host galaxies. An important, unresolved 
issue is whether age differences exist between the GC subpopulations in galaxies, and how large 
such age differences might be. This question is still controversial even in our Galaxy, where GC 
ages can be derived from color-magnitude diagrams (Stetson, Vandenberg, & Bolte 1996; Buonanno 
et al. 1998; Rosenberg et al. 1999), so it is not surprising that it has proven to be very difficult 
to obtain reliable information about ages of extragalactic GCs from integrated photometry and/or 
spectroscopy. As shown by de Freitas Pacheco & Barbuy (1995), blue horizontal-branch morphol- 
ogy can have a significant impact on observed H/3 indices, the most commonly used age indicator. 
More recently, Schiavon et al. (2002) have discussed a number of problems in the modeling of glob- 
ular cluster spectra, and concluded that uncertainties in the luminosity function of the red and 
asymptotic giant branches, a-element enhancement, horizontal branch morphology and other fac- 
tors can affect spectroscopic age estimates by at least 3-4 Gyr. Observationally, useful constraints 
on typical GC ages require very accurate measurements of Balmer line indices, as the separation 
between e.g. the 10 Gyr and 11 Gyr isochrones is only about 0.1 A in the equivalent width of H/?. 
Furthermore, at lower metallicities ([Fe/H] < —1), recent models (Maraston 2001, priv. comm., 
Bruzual & Chariot 2001, priv. comm.) suggest that the Balmer lines no longer depend on age in a 
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simple way for ages greater than about 10 Gyr. 

A more indirect way to constrain the formation history of stellar populations in general, in- 
cluding GCs, is by studying their abundance patterns. The Milky Way bulge and halo populations, 
as well as many early-type galaxies, show an enhancement in the abundances of a-elements relative 
to the Fe-group compared to the Sun (van den Bergh 1995, and references therein). The enhance- 
ment is typically around [a/Fe] ~ +0.3 (Carney 1996) and is also seen in the Sombrero bulge stars 
(Peletier et al. 1999). The most commonly adopted explanation is that different elements are pro- 
duced by stars with different lifetimes, of which supernovae of type la are generally assumed to have 
the longest lifetimes and make the strongest contribution to Fe. Thus, populations that formed on 
time-scales too short for SN la to have made any significant contribution to the abundance pat- 
terns show "enhancement" of other elements relative to Fe (e.g. Tinsley 1979). Therefore, studies of 
abundance patterns may provide an alternative chronometer for the chemical evolutionary history 
of galaxies, although the timescale for development of SN la remains very uncertain. 

Because of its relatively small discance and rich GC system, the Sombrero galaxy is an at- 
tractive target for detailed studies of individual GCs. With 8-10 m class telescopes, high-quality 
spectra can be obtained for many of the brighter clusters, and by co-adding spectra of several GCs 
it is possible to achieve sufficient S/N to bring the errors down to a level where model uncertainties 
rather than observational errors dominate. This makes it possible to test whether the models yield 
internally consistent results, e.g. by comparing age estimates based on different Balmer lines (H/3, 
H 7 , US). 

Spectroscopically, the GC system of the Sombrero has been studied by Bridges et al. (1997) 
who obtained rather low S/N spectra for 34 confirmed GCs. By co-adding the spectra they found 
a mean metallicity of [Fe/H] = —0.70 ± 0.3, but the S/N was too low to constrain metallicities of 
individual clusters. In this paper we present new spectra for GCs in the Sombrero, selected from 
the Paper I sample. Of 17 GC candidates, 14 are confirmed as GCs. Combining our data with 
those of Bridges et al. (1997), the total number of confirmed GCs in the Sombrero is increased to 
48. The best of our spectra have a S/N of about 50 per pixel, allowing us to accurately measure the 
metallicities of individual clusters and compare with previous photometric estimates. By co- adding 
the spectra of blue (metal-poor) and red (metal-rich) clusters, we also put constraints on the ages 
(Section 4.2) and abundance ratios (Section 4.3) of the Sombrero GC subpopulations. Finally, we 
use all available GC radial velocities in the Sombrero to constrain the mass of the Sombrero galaxy 
(Section 5). 



2. Data 

The data were obtained with the LRIS spectrograph (Oke et al. 1995) on the Keck I telescope 
during two observing runs on 2001, April 28 and May 19-20. We utilized the newly installed blue 
side on LRIS ( "LRIS-B" ) and a dichroic splitting at 5600A to simultaneously obtain blue and red 
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spectra, covering the regions ~3500-5500A and ~5700-8500A, respectively. On the blue side a 400 
1/mm grism was used, providing a dispersion of 1.74A pix -1 and a resolution of ~ 8A. On the red 
side we used a 600 1/mm grating blazed at 5000 A, with a dispersion of 1.28 A pix -1 and a resolution 
of ~ 6A. In addition to the multislit spectra of Sombrero GC candidates, we also obtained long-slit 
spectra for a number of radial velocity standard stars, flux standards, and Lick/IDS standard stars 
from Worthey et al. (1994). 

Because of anticipated problems with "ghosts" in LRIS-B, we applied 2" offsets to the telescope 
pointing for some exposures. Although these offsets were intended to be parallel to the slitlets, we 
found that they inevitably included some component perpendicular to the slits as well, resulting 
in greatly reduced signal from our science targets. Thus, we quickly decided to abandon this 
procedure. Another problem was related to the input coordinates. For the April run we used 
coordinates measured directly on the "central" and "inner halo" WFPC2 images from Paper I, 
using the METRIC task in the STSDAS package. Unfortunately, there were systematic differences 
between the two WFPC2 pointings which made it impossible to align all objects in the slitlets at 
the same time. For the May run we used direct images from the April run to alleviate this problem, 
but as a result of these slitmask alignment problems, we achieved a significantly lower S/N for some 
of the objects than could have been possible under optimal circumstances. 

Following initial processing of the images (bias subtraction, natfielding) , the spectra were 
extracted using standard tools in the IRAF SPECRED package 2 . As the spectra were obtained 
over several hours, significant shifts in the wavelength scale due to mechanical flexure were found to 
be present between individual integrations, amounting to up to ps 5 A. Assuming that all spectra in 
a given exposure were subject to the same shifts, the shifts were determined by cross-correlating the 
spectra of the brightest objects in each exposure with the spectra in a reference exposure, taken 
immediately before the arc spectra. For the red spectra, the wavelength zero-point was further 
fine-tuned using the NaD night sky lines at AA 5890 and 5896 A. Finally, the individual spectra of 
each object were co-added, using a sigma-clipping algorithm to reject cosmic ray events. The S/N 
at each pixel in the combined spectra was estimated from the variance of the combined pixels. The 
total integration time was 5^ hours. 

Table 1 lists all science objects for which spectra were obtained. The object IDs refer to Paper 
I, where the coordinates are listed. For convenience, we have also included the V, I photometry from 
Paper I. Note that the brightest cluster, H2-22, was saturated on the WFPC2 data used in Paper 
I, making its V —I color unreliable. The radial velocities in Table 1 are generally measured on the 
red spectra, for which the wavelength calibration zero-point could be checked using the night sky 
lines. The errors on the radial velocities in Table 1 are random errors, based on the uncertainties 
in the location of the cross-correlation peak reported by the FXCOR task. These do not take 
into account various systematic errors, mostly due to uncertainties in wavelength calibration. We 



2 IRAF is distributed by the National Optical Astronomical Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc. under contract with the National Science Foundation 
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estimate that such errors are on the order ~ 0.5A, or ~ 25 km/s. No skylines were present in 
the blue spectra, and differences of 100-150 km/s between radial velocities measured on the blue 
and red spectra were not uncommon. Another problem with the wavelength calibration in the 
blue stems from the fact that the calibration lamps in LRIS only had a few (5-7) suitable lines. 
The best fitting 2nd order polynomial wavelength solutions in the blue generally had a significant 
residual scatter (1-2A), presumably indicating that higher-order terms ought to be included in 
the wavelength solution, but the small number of calibration lines made it impossible to include 
higher-order terms. 

The radial velocity of the Sombrero galaxy is listed as 1024 km/s in the RC3 catalogue (de 
Vaucouleurs et al. 1991). If we exclude the three objects in Table 1 with RV<600 km/s, we are 
left with a sample of 14 globular clusters. These have a mean radial velocity of 1019 km/s and a 
dispersion of 246 km/s, in excellent agreement with the radial velocity of the Sombrero itself. Of the 
three objects with radial velocities less than 600 km/s, two (H2-06 and H2-10) are very compact 
and essentially unresolved on the HST images with effective radii < 0.3 pc (Larsen, Forbes & 
Brodie 2001) and are most likely foreground stars. The third (H2-27) is quite well-resolved on the 
HST images and has an effective radius of 3.4 pc, typical for globular clusters. It therefore seems 
likely that this object is indeed a globular cluster. The S/N of this spectrum is very poor and the 
cross-correlation signal was not very strong, so the radial velocity may be in error. In fact, using 
the blue spectrum of this object, we get a radial velocity of 1069 km/s. Because of the uncertain 
radial velocity and poor S/N we exclude this object from further analysis in this paper. 

Before spectrophotometric indices were measured, the spectra must be corrected for radial 
velocity. Rather than blindly applying the radial velocities in Table 1 to the blue spectra, we used 
a number of easily identifiable features such as the H/?, H7, US, Ca II H+K lines and the G-band 
to determine small zero-point shifts and scale corrections to the dispersion solutions. Thus, the 
wavelength scale of the final spectra used for the spectrophotometric analysis was typically accurate 
to better than lA in the region of interest (i.e. ~ 3700-5400 A). 

3. Measuring Lick/IDS indices on LRIS-B spectra 

In order to test if there might be any systematic differences between indices measured on LRIS- 
B spectra and the standard Lick/IDS system, we observed a number of stars from Worthey et al. 
(1994) in longslit mode. In Fig. 1 we compare our LRIS-B measurements of all Lick/IDS indices 
that were covered by our blue spectra with the standard values. Generally, our measurements agree 
well with the standard values, although there are slight systematic offsets for some of the narrower 
indices like Ca4227 and Ca4455. This may be attributed to the slightly higher resolution of LRIS-B 
compared to the Lick/IDS spectrograph. 

Unfortunately, the wavelength range covered in longslit mode by LRIS is slightly different from 
the one covered in multislit mode and does not extend beyond 5300A. Therefore, several key indices 
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(such as Mg2, Fe5270 and Fe5335) redward of this limit could not be compared. Furthermore, since 
this part of the spectrum is not covered by the flux standards, the flux calibration of the multislit 
spectra beyond 5300A also becomes uncertain. This is particularly problematic because the blue 
reflection efficiency of the dichroic already begins to drop around ~5300A. 

4. Constraining ages and metallicities 
4.1. Metallicities 

We derived metallicities for each GC using the calibration by Brodie and Huchra (1990), but 
adopting the modification to the error on the combined metallicity estimate described in Larsen 
and Brodie (2000). Metallicities for each of the individual calibrators, as well as the combined 
metallicity estimates are listed in Table 2. For comparison we have also included photometric 
metallicities based on the V—I colors in Table 1 and the calibration in Kissler-Patig et al. (1998). 
The spectroscopic and photometric metallicity estimates are compared in Fig. 2 and generally 
agree quite well. Three representative globular cluster spectra with low, intermediate and high 
metallicities are shown in Fig. 3. 

4.2. Comparison with population synthesis models 

High S/N spectra have the potential to remove the age-metallicity degeneracy inherent in 
photometric studies. At least in principle, ages are well constrained by Balmer line strengths, while 
metallicity information is provided by indices such as Fe5270, Fe5335 and Mg2 (Worthey 1994). 

Here we compare our Sombrero spectra with two sets of recent population synthesis models, 
kindly provided by C. Maraston (see Maraston & Thomas 2000) and R. Schiavon (Schiavon et al. 
2002b). Like most of the current models, the Maraston models are based on the "fitting" functions 
by Worthey et al. (1994). The Schiavon models are based on a new set of fitting functions derived 
from the spectral library of Jones (1999), and also include empirical corrections to the luminosity 
function of red giants (Schiavon et al. 2002). We use a preliminary version of the Schiavon models 
where only the Balmer line indices are based on the new fitting functions, i.e. other indices such 
as Mg6, Mg2 and the various Fe indices are still based on the Worthey fitting functions. The 
Schiavon models also cover a more limited age/metallicity range than the Maraston models, but the 
comparison nevertheless provides some useful insight into the differences between various available 
models. As discussed in Sect 3, our instrumental system appears to be very similar to the Lick/IDS 
system so we have not applied any corrections to spectrophotometric indices measured on the 
LRIS-B spectra. 

As mentioned earlier, some indices may be compromised by difficulties with the calibration 
of the blue spectra redward of 5300 A. The Mg2 index is particularly vulnerable to such problems 
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because of the wide separation between the continuum passbands, and offsets are frequently found 
between the standard Lick/IDS indices and instrumental systems (e.g. Kuntschner 2000). In Fig. 4 
we compare our measurements of Mg2 with Mgfr, which is slightly narrower but measures essentially 
the same feature. However, while the continuum passbands of Mg2 are located about 150A from the 
feature passband, those of Mgfe are located adjacent to it, making Mg6 much less sensitive to errors 
in the continuum slope. Generally, Fig. 4 shows that the two Mg sensitive indices correlate well, 
but the slope of our observed Mg6-Mg2 relation is shallower than the Maraston model predictions, 
which in turn agree well with Mg6-Mg2 plots in e.g. Trager et al. (1998) and Beasley et al. (2000). 
One possible explanation for this discrepancy is that errors in the continuum level at the position 
of the red continuum passband cause us to measure too low values for the Mg2 index because of 
uncorrected reduced throughput by the dichroic. Another possibility is that peculiarities in the 
abundance patterns affect the continuum passbands of Mg6 and Mg2 in different ways (Tripicco k, 
Bell 1995). 

In Fig. 5 we compare various index measurements with the Maraston and Schiavon model 
predictions. The Schiavon models are tabulated at 5, 7.9, 11.2 and 14.1 Gyr and at metallicities of 
[Fe/H] = —0.7, —0.4, 0.0 and +0.2. While the Maraston models are tabulated at 1 Gyr intervals 
from 2 to 15 Gyrs, we have plotted only the 5, 8, 11 and 14 Gyr models, for easy comparison 
with the Schiavon models. Note that the metallicity divisions of the Maraston models are different 
from those of the Schiavon models ([Fe/H]= -2.25, -1.35, -0.33, 0.00 and +0.35). Panels (a) and 
(b) show H/3 vs. (Fe) and the Mg6 index, where (Fe)= (Fe5270 + Fe5335)/2. The Maraston and 
Schiavon model grids are shown with solid lines and dashed lines, respectively. In panels (c) and 
(d) we show li^ and H5a vs. (Fe), using the definitions in Worthey & Ottaviani (1997). Because 
of the difficulties with the calibration of Mg2 discussed earlier, we have chosen Mgfo instead of 
Mg2 even though the formal errors on Mg6 are somewhat larger. The errors on the indices were 
estimated directly from the S /N of the spectra. 

To first order, Fig. 5 indicates high ages for most clusters, although the error bars are generally 
still too large to put strong constraints on the ages of individual clusters. In order to compare 
the mean properties of metal-rich and metal-poor GCs, we therefore co-added all spectra with 
[Fe/H]< —1 and [Fe/H]> —1, respectively. The co-added spectra are shown in Figure 6 and the 
corresponding datapoints are shown with filled circles in Fig. 5. In the H/3 plots both the metal- 
poor and metal-rich co-added datapoints fall near the ~ 11 Gyr Maraston isochrones, with little 
evidence of any significant age differences between the two within the formal ~ 2 Gyr error bars 
on the co-added data. The Schiavon models tend to predict lower ages for a given H/3 index than 
the Maraston models, especially at younger ages and lower metallicities. At the location of the 
metal-rich Sombrero GCs in the (Fe)- H/3 plane the difference is relatively modest and amounts to 
an age difference of only 1-2 Gyrs, but at lower metallicities the difference between the two sets of 
models becomes more pronounced and amounts to about 4 Gyrs at the lowest metallicity tabulated 
by Schiavon. Thus, the trend in these models is to decrease the ages of the metal-poor clusters 
more than those of the metal-rich ones, relative to the Maraston models. 
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The H7A and R5a measurements basically confirm the ages estimated from H/3, but note that 
both Hja and H5a reach a minimum around 11 Gyr at low metallicities and then increase strongly 
at higher ages. At [Fe/H] = —1.35, the Maraston models predict about the same EFy^ and H5a 
strengths at 15 and 5 Gyrs. The Schiavon models do not include tabulations of £[7^4 and H5a, 
but in Figure 7 we show their predictions for another H7 index defined by Vazdekis et al. (2001). 
This H7 index (in the following denoted H7 vaz ) differs from Hja in that the continuum passbands 
are located closer to the feature passband, thereby avoiding the strong metallicity dependence of 
H7A which results from its blue continuum passband coinciding with the G-band. In fact, Figure 7 
shows that H7 vaz is practically independent of metallicity, at least within the range spanned by 
the models. Taken at face value, Figure 7 would suggest that the metal-rich clusters are slightly 
older than the metal-poor ones, but both sub-populations are compatible with the same (high) age 
within the errors. Thus, although it should be stressed that even relative ages may still be somewhat 
model-dependent, it seems reasonable to conclude that both GC sub-populations in the Sombrero 
are old, and probably coeval to within the measurement uncertainties, i.e. with age differences of 
2-3 Gyrs at most. 



4.3. Abundance ratios 

Some attempts to model the effect of a-element enhancement on Lick indices have been made 
by Milone, Barbuy, & Schiavon (2000), who modeled the effect of a-element enhancement on the 
Lick Mg2 index and a broad index (TiOi2.5) covering the TiO band centered at 7257 A. The TiOi2.5 
index is included in our red LRIS spectra. Because of the problems with the calibration of Mg2, 
we converted the Milone et al. model predictions for Mg2 to Mgb, using the Maraston models. We 
found that the linear relation Mg6 = 14.03 Mg2 + 0.361 approximated the Maraston models to 
within 0.1 A for the age range 10-15 Gyr and Mg2 < 0.3 (Fig. 4). 

In addition to the strong terrestrial O2 absorption bands, which are located just outside the 
TiOi2.5 continuum passbands, the wavelength region around 7200-7500A is affected by terrestrial 
H2O absorption. We corrected our red spectra for this by normalizing the flux standard spectra with 
a smooth continuum and then dividing the GC spectra by the normalized flux standard spectra. 
This correction affected the TiOi2.5 index downwards by about 7A. Fig. 8 shows our measurements 
of TiOi2.5 an< i Mg6 on the co-added spectra, compared to the alpha-enhanced models by Milone, 
Barbuy, & Schiavon (2000). In the (Fe) — Mgb plot we also show the Maraston models for an 
age of 12 Gyr. Note that the most metal-poor Maraston models are actually rather similar to 
the [a/Fe] = +0.3 models by Milone, Barbuy, & Schiavon (2000), but approach the [a/Fe] = 
models at higher metallicities. This is probably because many of the metal-poor stars used to 
set up the Worthey fitting functions are alpha-enhanced (Maraston & Thomas 2000). For the 
metal-rich Sombrero GCs, both Mgfe and TiOi2.5 indicate clearly supersolar [a/Fe] ratios around 
[a/Fe] ~ +0.4. For the metal-poor clusters, the TiOi2.5 index again indicates supersolar [a/Fe] 
ratios similar to those of the metal-rich clusters, although the Mgb plot actually suggests [a/Fe] 
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closer to for the metal-poor clusters. However, at low metallicities the separation between the 
different [a/Fe] models is smaller and it seems most likely that both the metal-poor and metal-rich 
GCs in the Sombrero have super-solar [a/Fe] ratios similar to those observed in Milky Way GCs, 
the bulge of the Sombrero and in early-type galaxies. 



5. Mass of the Sombrero 

Globular clusters can be used as tracers of the mass of their host galaxy, as demonstrated e.g. 
by Huchra & Brodie (1987) (for M87) and Kissler-Patig et al. (1998) (NGC 1399). We combined 
our data in Table 1 with data for 34 Sombrero GCs from Bridges et al. (1997). Our sample does 
not overlap with that of Bridges et al. (1997), so by combining theirs and our data we get a total 
of 48 GCs. 

Fig. 9 shows the radial velocities as a function of the offset in R.A. from the center of the 
Sombrero (the major axis is aligned roughly east-west). Bridges et al. (1997) found hints of rotation 
in the cluster system, but our data clearly do not add much to the significance of their result. For 
completeness we also plot the radial velocities for the metal-rich and metal-poor clusters from our 
sample separately, but the data are clearly insufficient to provide any constraints on rotation in 
either sample. 

Following Huchra Sz Brodie (1987), we can use the GC radial velocities to obtain a mass 
estimate for the Sombrero in two ways. The virial mass is estimated as 

M VT = S ^ (1) 

and the projected mass is 

In these equations, is the distance between the ith and j'th clusters, is the distance between 
the ith cluster and the galaxy center and Vi is the velocity difference between the ith cluster and 
the mean system velocity. N is the number of clusters. The value of f p depends on the cluster 
orbits; following Bridges et al. (1997) we adopt f p = 16/ir for isotropic orbits. 

For the total sample of 48 clusters we get M V t = (8.1 ± 1.8) x 10 11 M and M P = (5.3 ± 
1.0) x 10 11 M within 17 kpc, adopting a distance of 8.7 Mpc (Larsen, Forbes & Brodie 2001). 
Using only our 14 clusters, we get M VT = (3.2 ± 1.5) x 10 11 M and M P = (4.2 ± 2.1) x 10 11 M 
(within 13 kpc). For comparison, Bridges et al. (1997) got M P = 5.2|;| x lO n M . The 12 clusters 
in the central HST pointing are useful as probes of the mass located within the central regions of 
the galaxy, and yield M VT = (2.1 ± 1.1) x 10 11 M and M P = (2.0 ± 1.3) x 10 11 M within the 
central 4.5 kpc. The errors were estimated by a standard "jackknife" procedure. 
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6. Summary and conclusions 

We have analyzed Keck / LRIS-B spectra for 17 globular cluster (GC) candidates in the 
Sombrero galaxy, of which 14 are confirmed as true GCs based on their radial velocities. One 
additional object has too low S/N for reliable radial velocity measurements, but based on the HST 
images presented by Larsen, Forbes & Brodie (2001) it is probably also a GC. We find generally 
good agreement between spectroscopic metallicities and previous photometric metallicites based on 
V — I colors. Comparison of Lick/IDS indices with population synthesis models indicates that both 
the metal-rich and metal-poor clusters are 10-15 Gyr old, with preferred ages around 11-12 Gyr. 
We see no detectable age differences between the two GC populations, although differences of a 
few Gyrs cannot be ruled out. Comparison of (Fe), Mg6 and TiO indicates that the Sombrero GCs 
probably have super-solar [a/Fe] ratios ~ +0.4, similar to or perhaps slightly larger than those 
observed in Milky Way globular clusters. 

Overall, the Sombrero GC system shares many similarities with that of the Milky Way, as 
well as many early-type galaxies. With few exceptions, both GC populations generally appear 
to be old, and another common feature seems to be an enhancement in a-element abundances, 
probably indicating that the gas out of which the clusters formed was enriched predominantly 
by type-II supernovae. This suggests that the majority of GCs formed early on, before type la 
SNe made a significant contribution to chemical enrichment, and presumably within a few Gyrs 
after star formation commenced. The metal-rich GCs in the Sombrero, like those in most galaxies, 
probably formed early on when galaxies were still assembling with high star formation rates and gas 
densities providing favorable conditions for the formation of large numbers of massive star clusters. 
The metal-poor clusters might have formed even earlier at high redshifts, triggered by ionization 
fronts (Cen 2001). 

Finally, we combine our sample with data for 34 GCs from Bridges et al. (1997) and obtain 
virial and projected masses for the Sombrero of Mvt = (8.1 ± 1.8) x 10 11 M and Mp = (5.3 ± 
1.0) x 1O 11 M , respectively. This is in good agreement with previous estimates of the mass of 
the Sombrero. For 12 clusters located within the central 4.5 kpc (projected) we estimate Mvt = 
(2.1 ± 1.1) x 10 11 M and M P = (2.0 ± 1.3) x 10 11 M . 

This work was supported by National Science Foundation grant number AST9900732. We 
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Dr. Ricardo Schiavon for allowing us to use his population models prior to publication. Helpful 
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Fig. 1. — Lick/IDS indices measured on our LRIS-B spectra, compared to standard values from 
Worthey et al. (1994). The agreement is generally good, except for narrow indices like Ca4227 and 
Ca4455. 
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Fig. 2.- 



Comparison of spectroscopic and photometric metallicity estimates for Sombrero GCs. 
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Fig. 3. — Spectra for three globular clusters with low (top), moderate (center) and high (bottom) 
metallicities. Note: no smoothing has been done. The two upper spectra have been shifted for 
clarity, with the zero- levels indicated by the short dashed lines. 
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Fig. 4. — Comparison of Mg6 and Mg2 for Sombrero GCs. Maraston models for 10 and 15 Gyr arc 
overplotted on the figure. 
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Fig. 5. — Comparison of our data with population synthesis models by C. Maraston and R. Schi- 
avon. Panels (a) and (b): (Fe) and Mgb vs. H/3. Panels (c) and (d): (Fe) vs. H^a an d HSa- 
Maraston models (solid line grids) are shown for ages of 5, 8, 11 and 14 Gyr and [Fe/H] = —2.25, 
— 1.35, —0.33, 0.0 and +0.35. Schiavon models (dashed line grids) are for ages of 5, 7.9, 11.2 
and 14.1 Gyr and [Fe/H] = —0.7, —0.4, 0.0 and +0.2. The filled circles represent the co-added 
metal-poor and metal-rich spectra. 
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Fig. 6. — Co-added spectra of all metal-rich ([Fe/H] > —1.0) and metal-poor ([Fe/H] < —1.0) 
globular clusters. Panels (a) and (b) show the spectra from the blue and red side arms on LRIS, 
respectively. Lick indices used in this paper are indicated. The metal-poor spectra have been 
shifted upwards, as indicated by the horizontal dashed lines. 
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Fig. 7. — Measurements of H7 vaz according to the Vazdekis et al. (2001) definition and (Fe), 
compared to Schiavon (2002) models. 
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Fig. 8. — TiOi2.5 an d Mgb vs. (Fe) for co-added metal-rich and metal-poor spectra, compared to 
a-enhanced models by Milone, Barbuy, & Schiavon (2000). The dashed line in the Mg6 vs. (Fe) 
plot represents 12 Gyr Maraston models. 



- 21 - 



600 
400 
200 


-200 
-400 



600 - 
400 r- 
200 r- 

o '-■ 

-200 r- 
-400 r- 
600 F~ 
400 r- 
200 r- 


-200 r- 
-400 '- 



All 



■<>■ 



*1. % 



31 - 

2fc 



. 2jS — 

as 



Metal -rich 



Metal-poor 



-4 -2 2 4 

R.A. offset (orcmin) 

Fig. 9. — Radial velocity as a function of distance from the center along the major axis for our 
full GC sample combined with data from Bridges et al. (1997) (top), our metal-rich clusters only 
(center) and our metal-poor clusters only (bottom). Data from Bridges et al. (1997) are shown 
with open circles while our data are shown with filled circles. No significant evidence for rotation 
is seen in either sample. 
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Table 1. Globular cluster candidates in NGC 4594 



ID 


V 


fir 




KV 

km/s 


Q /AT 
b/JN 

Blue F 


led 




1 / H2-06 


20.28 





.92 


304 ± ll b 


46 




N 


2 / H2-10 


20.51 


1 


.01 


108 ±8 


45 


55 


N 


4 / H2-09 


20.45 





.97 


611 ±9 


45 


44 


Y 


5 / H2-22 


19.1 a 






1274 ± 17 


19 


24 


Y 


6 / H2-27 


22.04 


1 


.14 


227 ±39 


9 


9 


Y? 


7 / C-076 


21.36 





.90 


1563 ± 77 b 


8 


16 


Y 


8 / C-068 


20.37 


1 


.11 


849 ± 11 


32 


41 


Y 


9 / C-051 


21.17 


1 


.18 


1015 ± 18 


17 


25 


Y 


10 / C-064 


21.10 





.94 


865 ± 16 


17 


23 


Y 


11 / C-059 


21.14 


1 


.17 


1034 ± 19 


13 


22 


Y 


12 / C-032 


20.38 





.97 


1055 ± 24 


25 


36 


Y 


13 / C-042 


21.19 


1 


.26 


1225 ±26 


11 


15 


Y 


15 / C-116 


19.98 





.96 


1022 ± 8 


46 


46 


Y 


16 / C-136 


21.53 


1 


.06 


1202 ± 17 


14 


14 


Y 


17 / C-132 


20.91 


1 


.15 


679 ± 15 


16 


23 


Y 


18 / C-134 


20.29 


1 


.24 


950 ± 15 


41 


34 


Y 


19 / C-137 


21.10 


1 


.19 


919 ±30 


24 


33 


Y 



Note. — RV = heliocentric radial velocity, V = visual magnitude. The S/N column gives the average 
signal-to-noise per pixel in the region 4500A- 5000A (blue) and 6500A- 7000A(red). a : This object was 
saturated on the WFPC2 data used in Paper I. b : Radial vel. from blue spectrum. The IDs refer to Larsen, 
Forbes & Brodie (2001), where coordinates are listed. The quoted errors on the radial velocities are random 
errors determined from the uncertainty on the cross-correlation peaks as reported by FXCOR. Additional 
systematic errors are ~ 25 km/s (see text). 



Table 2. Metallicity estimates based on the calibration in Brodie and Huchra (1990). The two 
bottom lines give data for the co-added metal-poor ([Fe/H] < —1) and metal-rich ([Fe/H] > — 1) 

spectra. 



Object 








Individual [Fe/H] estimators 








[Fc/H] (sp) [Fc/H] (phot) 




A 


Mg 2 


MgH 


Gband CNB 


Fe5270 


CNR. 


H+K 





4 


-1 


.59 


± 





37 


-2 


2!) 


± 





34 


-2 


.35 


± 





.48 


-1. 


72 


± 





.33 


-1 


.96 


± 





34 


-1 


57 


± 





62 


-1 


45 


± 





45 


-1 


.72 


± 





43 


-1 


86 


± 


0. 


14 


-1. 


23 


5 


-1 


28 


± 





.37 


-2 


23 


± 





.36 


-2 


.06 


± 





.52 


-1 


07 


± 





41 


-1 


.34 


± 





41 


-1 


63 


± 





.65 


-1 


.30 


± 





49 


-0 


19 


± 





51 


-1 


41 


± 


0. 


27 


-0. 


42 


7 


-0 


.87 


± 





.39 


-2 


10 


± 





.43 


-1 


.30 


± 





.67 


-1. 


85 


± 





.64 


-0 


.78 


± 





.97 


-1 


98 


± 





.82 


-2 


40 


± 





54 


-2 


92 


± 


1 


.28 


-1 


.80 


± 


0. 


31 


-1. 


50 


8 


-0 


47 


± 





37 


-1 


10 


± 





.34 


-1 


17 


± 





.49 


-0. 


50 


± 





.35 


-0 


.97 


± 





36 


-0 


90 


± 





63 


-0 


96 


± 





46 


-0 


26 


± 





45 


-0 


81 


± 


0. 


14 


-0. 


67 


9 





.04 


± 





37 


-0 


.68 


± 





.37 


-0 


.73 


± 





.53 


-1. 


80 







43 





.01 


± 





48 


-0. 


60 


± 





70 


-0 


.63 


± 





48 





55 


± 





.53 


-0 


54 


± 


0. 


31 


-0. 


52 


10 


-1 


.06 


± 





.37 


-2 


11 


± 





.36 


-1 


22 


± 





.54 


-1 


33 


± 





.40 


-1 


.17 


± 





.35 


-1 


60 


± 





.68 


-1 


42 


± 





47 


-1 


44 


± 





48 


-1 


44 


± 


0. 


15 


-1 


21 


11 





24 


± 





.37 


-0 


84 


± 





.40 


-1 


84 


± 





.59 


-0. 


86 


± 





.48 


-0 


.63 


± 





59 


-0 


58 


± 





.76 


-0 


.56 


± 





47 


-0 


56 


± 





61 


-0 


74 


± 


0. 


23 


-0. 


62 


12 


-0 


.90 


± 





37 


-1 


57 


± 





35 


-1 


24 


± 





.51 


-0. 


45 


± 





.35 


-1 


.42 


± 





.35 


-1 


58 


± 





.65 


-1 


24 


± 





46 


-0 


.38 


± 





.45 


-1 


.09 


± 


0. 


20 


-1 


23 


13 


-0 


.59 


± 





.38 


-1 


.03 


± 





41 


-1 


91 


± 





.60 


0. 


16 


± 





.68 


-1 


.05 


± 





44 


-1 


04 


± 





83 


-0 


.73 


± 





.53 


-0 


.59 


± 





66 


-0 


.88 


± 


0. 


24 


-0. 


23 


15 


-0 


.88 


± 





.37 


-1 


62 


± 





.34 


-1 


.83 


± 





.49 


-1 


08 


± 





.32 


-1 


.02 


± 





34 


-1 


17 


± 





.62 


-1 


26 


± 





45 


-0 


.59 


± 





43 


-1 


19 


± 


0. 


16 


-1 


22 


16 


-1 


.08 


± 





37 


-1 


62 


± 





.39 


-2 


34 


± 





.58 


-1 


09 


± 





.46 


-1 


.43 







39 


-1 


56 


± 





81 


-1 


.32 


± 





48 


-0 


84 


± 





.53 


-1 


41 


± 


0. 


18 


-1 


10 


17 


-0 


.30 


± 





37 


-0 


.32 


± 





.38 


-1 


79 


± 





.55 


-0. 


52 


± 





48 


-0 


.75 


± 





42 


-1 


02 


± 





71 





.20 


± 





49 





36 


± 





.53 


-0 


49 


± 


0. 


28 


-0. 


67 


18 


-0 


20 


± 





37 


-0 


.08 


± 





.35 


-0 


76 


± 





.49 


-0. 


24 


± 





.33 


-0 


.17 


± 





34 


0. 


33 


± 





.63 





42 


± 





46 


-0 


15 


± 





43 


-0 


12 


± 


0. 


14 


-0 


39 


19 


-0 


07 


± 





.37 


-0 


78 


± 





.36 


-1 


.00 


± 





.51 


-0. 


71 


± 





.46 





24 


± 





37 


-0 


90 


± 





66 


-0 


47 


± 





47 





11 


± 





47 


-0 


43 


± 


0. 


21 


-0. 


47 


metal-poor 


-1 


11 


± 





37 


-1 


82 


± 





.34 


-1 


81 


± 





.48 


-1 


16 


± 





.33 


-1 


.42 


± 





34 


-1 


44 


± 





61 


-1 


32 


± 





45 


-0 


90 


± 





43 


-1 


39 


± 


0. 


14 






metal-rich 


-0 


24 


± 





37 


-0 


67 


± 





.34 


-1 


17 


± 





.49 


-0. 


58 


± 





.35 


-0 


.54 


± 





35 


-0 


54 


± 





.62 


-0 


.35 


± 





46 


-0 


11 


± 





45 


-0 


54 


± 


0. 


13 







Table 3. Indices in the Lick/IDS system 



Object 


H/3 




H6 A 


H7 vaz 


Fe5270 


Fo5335 


(Fe) 


Mg 2 


Mgb 


MgFe 


Ti0 12 . 5 




A 


A 


A 


A 


A 


A 


A 


mag 


A 


A 


A 



4 


1 


74 


± 





.20 


1 


.59 


± 





.26 


2 


.25 


± 





.29 


1 


.34 


± 





24 





92 


± 





20 


1 


.15 


± 





.26 


1 


.03 







16 


-0 


.00 


± 


0. 


00 


1 


.01 


± 





.18 


1 


02 


± 


0. 


.12 


2.1 ± 7. 


4 


5 


1 


.69 


± 





.46 





.55 


± 





.74 


1 


04 


± 





.83 


1 


.03 


± 





.81 


1. 


.11 


± 





.42 


1 


.16 


± 





.53 


1 


.13 


± 





.34 





.01 


± 


0. 


01 


0. 


42 


± 





.45 


0. 


09 


± 


0. 


39 


-24.2 ± 13. 


6 


7 


1 


.86 


± 


1 


.17 


3 


.16 


± 


1 


.22 


5. 


11 


± 


1 


.60 


1 


.40 


± 


1 


58 


0. 


25 


± 





.97 





.05 


± 


1 


.36 





15 


± 





.84 





.03 


± 


0. 


03 


-0. 


25 


± 


1 


.11 


0. 


00 


± 


0. 


84 


-13.0 ± 15. 


5 


8 


1 


.61 


± 





.33 


-4 


.07 


± 





.37 


-0 


39 


± 





.34 


1 


.10 


± 





.29 


1 


89 


± 





.27 


2 


.73 


± 





.35 


2 


.31 


± 





.22 





.12 


± 


0. 


01 


2. 


95 


± 





.25 


2. 


.61 


± 


0. 


.17 


21.9 ± 6. 


7 


9 


1 


.95 


± 





.53 


-4 


12 


± 





.67 


-1. 


65 


± 





.86 





.17 


± 





.64 


2. 


93 


± 





.59 


2 


.94 


± 





.60 


2 


.94 


± 





.42 





17 


± 


0. 


01 


3. 


00 


± 







3. 


25 


± 


0. 


34 


39.4 ± 13. 


8 


10 


3 


.73 


± 





.53 


-0 


.20 


± 





.59 


2. 


20 


± 





.69 


1 


.71 


± 





.61 


1 


11 


± 





.54 


-0 


.68 


± 





.66 





.22 


± 





.42 





.03 


± 


0. 


01 


1 


20 


± 





.50 


0. 


51 


± 


0. 


51 


10.0 ± 16. 


3 


11 


2 


.59 


± 





.67 


-6 


.95 


± 





.86 


-3 


49 


± 





.35 





.37 ± 





76 


3. 


01 


± 





.77 





.48 


± 


1 


.11 


1 


.74 


± 





.68 





14 


± 


0. 


02 


4. 


00 


± 





.58 


2. 


04 


± 


0. 


55 


16.7 ± 15. 


1 


12 


2 


.85 


± 





.38 


-3 


04 


± 





.40 


1. 


00 


± 





42 





.75 


± 





42 


0. 


75 


± 





.39 


1 


.57 


± 





.52 


1 


.16 


± 





.33 





.08 


± 


0. 


01 


1 


24 


± 





.51 


1 


20 


± 


0. 


30 


8.6 ± 7. 


8 


13 


2 


.91 


± 





.82 


-3 


.85 


± 


1 


.27 


1 


65 


± 


1 


.46 





.85 


± 


1 


.31 


1. 


.05 


± 





.97 


1 


.90 


± 


1 


.16 


1 


.47 


± 





.75 





14 


± 


0. 


02 


2. 


09 


± 





.90 


1 


99 


± 


0. 


.61 


30.0 ± 24. 


6 


15 


2 


.10 


± 





.22 


-0 


.54 


± 





.21 


2. 


35 


± 





.21 


1 


.57 


± 





.18 


1 


70 


± 





.21 


1 


.84 


± 





.30 


1 


.77 


± 





.18 





.06 


± 


0. 


01 


1 


52 


± 





.22 


1 


.64 


± 


0. 


15 


28.3 ± 7. 





16 


3 


31 


± 





.65 


-1 


.33 


± 





.90 


2. 


83 


± 





.86 





.87 


± 





84 


0. 


90 


± 





.96 


1 


.75 


± 


1 


.25 


1 


33 


± 





.79 





.05 


± 


0. 


02 


2. 


11 


± 





.67 


1 


.67 


± 


0. 


57 


-6.4 ± 24. 


1 


17 


1 


61 


± 





.58 


-5 


.00 


± 





.66 


-3 


44 


± 





.94 


1 


.87 


± 





.02 


2. 


48 


± 





.67 





.87 


± 





.87 


1 


.68 


± 





.55 





.20 


± 


0. 
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